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PREFACE 


The  investigation  described  in  this  report  was  authorized  as  part  of  the 
Civil  Works  Research  and  Development  Program  by  the  Headquarters,  US  Army 
Corps  cf  Engineers  (HQUSACE) .  This  study  was  conducted  by  the  work  unit  enti¬ 
tled  "Stochastic  Coastal  Models,"  Work  Unit  32487,  under  the  Shore  Protection 
and  Restoration  Program  at  the  Coastal  Engineering  Research  Center  (CERC)  of 
the  US  Army  Engineer  Waterways  Experiment  Station  (WES).  Messrs.  John  H. 
Lockhart,  Jr.,  John  G.  Housley,  and  James  E.  Crews  were  the  HQUSACE  Technical 
Monitors.  Dr.  C.  Linwood  Vincent  is  CERC  Program  Manager. 

The  study  was  conducted  from  July  1986  through  November  1987  by 
Dr.  Todd  L.  Walton,  Jr.,  CERC,  under  the  general  supervision  of  Dr.  James  R. 
Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Chief  and  Assistant  Chief,  CERC, 
respectively,  and  under  the  direct  supervision  of  Mr.  Thomas  W.  Richardson, 
Chief,  Engineering  Development  Division,  CERC. 

Commander  and  Director  of  WES  during  publication  of  this  report  was 
COL  Larry  B.  Fulton,  EN.  Dr.  Robert  W.  Whalin  was  Technical  Director. 
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INTRODUCTION 


Considerable  debate  (CHRYSTALand  MURRAY, 1906;  MCDONALD, 1954;  BRUNK.1959 
PEC'COK,19«l;  MEGERIAN.1964;  LIU, 1970;  COHN  and  ROBINSON, 1976;  QUINN, 1986; 
CHANGNON.1987)  has  been  ongoing  since  the  early  portion  of  the  twentieth  cen¬ 
tury  concerning  man’s  ability  to  predict  potential  future  periods  of  time  during 
which  lake  levels  would  be  at  high  stages.  Storms  that  occur  during  these  periods  of 
high  lake  level  do  considerable  damage  to  shore  structures,  beaches,  bluffs,  and  lake- 
front  dwellings.  Attempts  to  mitigate  erosion  and  potential  flood  damage  from  high 
lake  stages  while  at  the  same  time  improving  navigation  benefits  by  prevention  of 
excessive  periods  of  low  water  levels  have  resulted  in  a  series  of  complex  water  level 
control  regulations  on  both  Lake  Superior  and  Lake  Ontario  (INTERNATIONAL 
JOINT  COMMISSION, 1976)  altering  natural  stages  of  the  lakes  to  some  extent. 
Additional  effects  of  channel  dredging  between  the  lakes,  as  well  as  flow  diversions 
to  and  from  the  lakes  have  also  altered  natural  lake  stages,  but  to  a  lesser  degree. 

The  extent  of  man’s  influence  on  the  lakes  is  subject  to  wide  debate  but  appears 
minor  compared  to  nature’s  cycles  or  pseudocyles  in  lake  water  level.  This  hypoth¬ 
esis  will  be  investigated  in  a  later  portion  of  this  paper  for  the  more  recent  period 
of  time  from  1921  thru  1987. 

An  ability  to  forecast  future  water  levels  would  be  useful  for  optimal  control  of 
lake  stages  but  it  is  also  a  necessary  first  step  in  simulating  future  water  levels  for 
prediction  of  potential  beach  erosion  that  might  occur  in  a  major  storm.  Operational 
forecast  models  of  lake  stages  appear  to  be  of  a  short  forecast  horizon  (DECOOKE 
and  MEGERIAN.1967;  CROWLEY, 1987)  with  the  apparent  belief  that,  the  natural 
driving  force  for  lake  level  change  being  random  precipitation,  long  range  (>  6 
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months)  forecasting  would  be  of  limited  utility.  It  is  important  to  note  though  that 
any  persistence  in  the  data  (i.e.  high  lake  levels  followed  by  high  lake  levels  or  vice 
versa)  justifies  an  attempt  to  extend  the  forecast  horizon  for  lake  stage  prediction  and 
hence  provide  better  ability  to  do  realistic  simulation  of  future  lake  level  events.  The 
present  paper  is  an  attempt  to  utilize  a  recent  geophysical  digital  signal  processing 
algorithm,  the  Maximum  Entropy  Method  (BURG,  1975)  to  forecast  Great  Lake 
monthly  average  water  levels  in  an  attempt  to  separate  monthly  average  water  level 
signal  informational  content  from  random  noise  in  the  time  series.  By  reducing  the 
water  level  time  series  to  random  uncorrelated  noise,  the  potential  for  simulation  of 
water  levels  containing  the  same  statistical  properties  and  correlation  structure  as 
the  original  series  is  realized.  The  methodology  used  to  do  the  forecast  is  shown  to 
provide  a  forecast  horizon  well  beyond  one  year  in  the  case  of  the  lower  Great  Lakes 
(Superior  is  intensively  regulated  and  therefore  lake  level  forecasts  are  very  close  to 
actual  with  little  noise  in  the  series). 

PAST  INVESTIGATIONS 

Various  investigators  have  attempted  to  predict  future  cycles  of  high  water  in 
the  Great  Lakes  for  the  previously  mentioned  reasons.  COHN  and  ROBINSON 
(1976)  provided  a  deterministic  process  Fourier  analysis  approach  to  the  problem  of 
determining  monthly  average  lake  levels  with  the  result  that  the  lake  levels  followed 
natural  cycles  with  periods  of  1,11,  22,  and  36  years.  L1U(1970),  in  a  more  recent 
analysis,  investigated  the  spectral  shape  of  the  monthly  average  water  levels  as  well 
as  the  probability  distributions  of  monthly  average  water  levels.  LIU(1970)  found 


that  the  water  levels  were  approximately  Gaussian  distributed,  but  often  failing  a 
Chi  Square  test  at  a  95  percent  significance  level.  He  noted  that  the  series  was  better 
fit  with  a  Gram-Charlier  probability  distribution  function  which  allows  for  skewness 
within  the  data,  and  he  attributed  the  skewness  to  possible  nonlinear  aspects  of  the 
hydrological  processes  involved  in  lake  level  forcing  functions.  As  the  Chi  squared 
test  is  a  weak  test  for  normality  (STEPHENS  (1974)),  the  hypothesis  of  normality 
will  later  be  tested  with  a  stronger  statistical  test,  the  Kolmogorov-Smirnov  test. 
LIU ’S(  1970)  conclusions  concerning  the  spectral  shape  of  the  deseasonalized  time 
series  suggested  that  the  lake  levels  were  indeed  noise  with  no  apparent  true  cycles 
except  for  the  seasonal  cycle.  An  important  point  not  mentioned  was  that  the  series 
was  correlated  noise,  not  white  noise,  hence  additional  informational  content  still 
existed  within  the  series  to  improve  upon  future  prediction  and  hence  simulation. 

YEVJEVICH  (1975), in  an  approach  similar  to  the  one  being  presented  here,  de¬ 
veloped  a  simulation  model  for  the  mean  monthly  net  basin  supply  for  the  various 
lakes,  where  the  net  basin  supply  consists  of  precipitation,  evaporation,  groundwater 
contribution,  surface  inflow  and  outflow  to  the  lake.  As  the  net  basin  supplies  were 
derived  indirectly  from  limited  noisy  measurements  and  incorporated  additional  un¬ 
certainty  in  estimating  evaporation  and  groundwater  contribution,  the  end  result 
was  a  noisy  signal  with  limited  informational  content. 

The  present  approach  to  future  water  level  simulation  was  to  utilize  the  exist¬ 
ing  data  base  on  Great  Lakes  water  levels  (NATIONAL  OCEANIC  and  ATMO¬ 
SPHERIC  ADMINISTRATION, 1987)  due  to  the  fact  the  gage  water  levels  present 
an  accurately  measured  time  series  relatively  uncontaminated  by  man’s  lack  of  un¬ 
derstanding  of  the  natural  physical  processes  of  precipitation  and  evaporation. 
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MODEL 

The  model  proposed  here  is  a  hybrid  model  composed  of  a  seasonal  component 
combined  with  a  linear  autoregressive  model,  the  Maximum  Entropy  Method  model 
(BURG, 1975).  This  model  will  be  shown  to  produce  a  reasonable  forecast  and  be 
capable  of  reducing  the  original  monthly  average  water  level  to  random  uncorrelated 
noise.  The  period  of  record  used  to  test  the  model  in  a  forecast  mode  was  1921  thru 
1987.  Although  Great  Lakes  monthly  average  water  levels  date  back  to  the  1860’s 
prior  to  regulation  of  Lake  Superior,  the  control  plan  and  operating  rules  initiated 
in  1921  on  Lake  Superior  appear  to  have  modified  the  lake  level  variations  on  the 
lower  Great  Lakes  (INTERNATIONAL  JOINT  COMMISSION, 1975).  Figures  1 
through  3  show  the  monthly  average  water  levels  for  gages  on  Lakes  Michigan-Huron 
(5014),  Lake  Erie  (3063),  and  for  Lake  Ontario  (2030)  for  the  period  1921  thru  1987. 
Autocorrelation  plots  for  the  same  series  ace  shown  in  Figures  4  through  6  and  show 
the  high  persistence  along  with  seasonal  correlation  within  the  lake  level  time  series. 
This  persistance  is  marked  by  not  only  a  strong  correlation  between  months  but  also 
a  strong  correlation  at  yearly  levels  (  lag  =  12  months).  This  particular  structure 
with  correlation  at  seasonal  (in  this  case  ye  My  lags)  is  suggestive  of  a  seasonal  trend 
(either  nonstationary  random,  stationary  random,  or  deterministic)  in  the  data  and 
can  be  handled  by  a  variety  of  methods. 

If  the  seasonal  trend  is  assumed  to  be  of  a  nonstationary  nature  (i.e.  does  not 
seek  a  seaspnal  time  dependent  level),  then  the  best  approach  to  deseasonalizing 
the  data  is  via  seasonal  differencing,  in  the  cue  of  monthly  average  data  this  in¬ 
volves  differencing  at  lag  12.  In  terms  of  a  linear  filter  representation,  this  type  of 
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deseasonalizing  filter  L  given  as: 

Z(<)  =  (1-S5)X(<)  =  V,X(<)  (1) 


where  X(t)  is  the  original  series,  Z(t)  is  the  deseasonalized  series  by  this  approach, 
and  B 3  wi:h  s=12  is  a  backwards  lag  12  operator  (i.e.  see  BOX  and  JENKINS 
(1976)).  The  term  1 1  —  B12)  =  V12  is  a  linear  difference  operator  or  filter  of  order 
1  and  degree  12  on  the  data.  BOX  and  JENKINS  (1976)  are  proponents  of  this 
approach  to  handling  seasonal  correlation  structure  in  the  data  and  recommend  sea¬ 
sonal  differencing  to  reduce  the  data  to  stationarity.  A  particular  problem  with  this 
approach  is  that  such  differencing  removes  a  stationary  mean  level  of  the  process 
and  hence  forces  a  nonstationarity  into  the  data  regardless  of  whet  he  the  data  are 
nonstationary  or  not.  This  type  of  filter  is  equivalent  to  an  assumed  random  walk 
(at  seasonal  lag)  model  assumption.  Although  nonstationarity  may  be  present  in 
long  geophysical  records  such  as  the  ones  analyzed  here  due  to  long  term  climatic 
changes  or  man  made  changes  (i.e.  flow  diversions,  channel  deepening,  or  urbaniza¬ 
tion  effecting  runoff  changes),  such  changes  should  not  be  assumed  a  priori  without 
physical  justification. 
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autoregressive  polynomial  are  all  less  than  one  (i.e.  inside  the  unit  circle  in  the 
compl 'X  plane)  is  equivalent  to  assuming  stationarity.  BOX  and  JENKINS  (1976) 
are  also  proponents  of  this  approach  to  modeling  and  recommend  deseasonalizing  of 
this  type  in  conjunction  with  similar  techniques  (i.e.  ARM  A  models)  for  analysis 
of  the  deseason  a  lized  series  When  applying  the  system  of  filters  as  given  in  Eq.  2 
along  with  a  stationary  data  filter  given  in  general  form  as: 


(3) 


M(B) 

where  W(t)  represents  a  wh:te  noise  (inno\ ation)  term  and  the  roots  of  the  autore¬ 
gressive  polynomial  are  less  than  absolute  value  one,  the  model  becomes  what  BOX 
and  J  ENKINS(1976)  refer  to  as  a  seasonal  multiplicative  model.  A  possible  criticism 
of  the  multiplicative  seasonal  model  is  that  there  is  no  physical  reason  to  believe  that 
a  white  noise  process  W(t)  should  be  influenced  by  antecedent  conditions  at  a  lag 
period  s.  This  particular  model  was  used  in  the  investigation  of  lake  levels  and  will 
be  discussed  further  below. 

A  third  approach  to  the  deseasonalization  of  the  level  data  is  provided  by  stan¬ 
dardization  of  the  seasonal  series  as  given  by: 

X  (<)  —  m(t) 


Z(t) 


Sd(t) 


(4) 


where  m(t.)  is  the  estimated  time  varying  mean  value  of  X(t),  in  the  present  case 
monthly  averages,  and  Sd(t)  is  the  time  varying  standard  deviation  of  X(t)  (i.e. 
monthly  indexed  standard  deviation).  This  particular  approach  is  well  adapted  to 
heteroscedastic  data  of  the  type  seen  in  hydrologic  time  series  as  noted  by  KASHYAP 
and  RAO(1976).  After  deseasonalization  by  this  approach,  the  data  typically  assume 
stationary  form  to  which  a  stationary  type  of  model  such  as  an  AR  or  ARMA 
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model  can  be  fit.  For  parsimony  of  terms,  the  time  varying  monthly  mean  and  vari¬ 
ance  terms  can  be  fit  by  discrete  Fourier  series  with  only  those  coefficients  retained 
which  show  significance.  In  the  present  analysis,  this  standardization  method  of 
dcscasonalization  was  used  but  the  seasonal  coefficients  were  used  (12  means  and  12 
standard  deviations)  rather  than  providing  a  Fourier  series  fit  to  the  seasonal  means 
and  variances. 

Three  methods  were  used  to  assess  whether  the  assumption  of  stationarity  was 
reasonable  for  the  water  level  monthly  average  series  during  the  period  of  record 
from  1921-1987.  The  first  test  used  was  a  runs  test  as  suggested  in  BENDAT  and 
PIF.RSOL(197l),  which  consists  of  counting  the  number  of  runs  about  the  mean  of 
the  process  for  the  mean  square  value  of  equal  segments  of  record.  This  particular 
test  does  not  require  that  the  seasonal  trend  be  removed  prior  to  testing  as  long  as 
the  fundamental  period  is  short  compared  to  the  averaging  time  used  to  compute 
sample  values,  therefore  the  series  was  not  deseasonalized  prior  to  testing.  The 
hypothesis  that  the  data  were  random  stationary  at  a  5  percent  significance  level 
was  tested  and  accepted  for  all  three  lake  level  series. 

Prior  to  the  remaining  two  methods  for  testing  stationarity,  the  data  were  desea¬ 
sonalized  via  standardization  with  monthly  means  and  variances  as  noted  in  previous 
paragraphs.  The  second  method  used  to  assess  stationarity  was  to  split  each  series 
into  two  segiments,  beginning  1921-1954  and  ending  1955-1987,  and  compare  both 
autocorrelation  functions  and  raw  spectrums  of  the  data  for  consistency  throughout 
time.  This  method  also  did  not  provide  any  evidence  for  nonstationarity  in  the 
series. 

The  third  method  utilized  for  checking  the  stationarity  assumption  was  tosegiment 
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the  series  as  in  the  method  above  and  compute  autoregressive  parameters  for  the 
series  utilizing  the  BURG(1967)  Maximum  Entropy  Method  algorithm.  Results  of 
this  particular  test  for  a  two  parameter  autoregressive  model  test  are  given  in  Table 
1.  Although  differences  in  the  parameter  estimates  between  the  first  half  1921-1954 
and  last  half  1955-1987  of  the  series  for  the  three  gages  tested  are  noted,  it  would 
appear  that  such  differences  are  of  minor  significance.  This  observation  is  based  on 
the  fact  that  the  first  autoregressive  parameter  estimate  varies  less  than  10  percent 
and  averages  only  5  percent  difference  between  the  two  record  halves. 

As  noted  in  the  previous  paragraph,  the  water  level  time  series  for  the  lower 
lakes  did  not  show  evidence  for  nonstationary  after  1921  (the  year  in  which  major 
regulation  of  Lake  Superior  was  initiated)  in  any  of  the  nonstationary  tests.  This 
analysis  suggests  that  man’s  intervention  on  the  lakes  has  not  seriously  altered  the 
natural  levels  of  the  lower  Great  Lakes  at  least  in  the  period  since  1921.  It  is  on 
the  basis  of  this  preliminary  data  checking  that  a  univariate  model  for  lake  level 
simulation  is  based.  As  the  series  appears  to  be  generally  stationary,  a  constant 
parameter  model  will  be  assumed  sufficient  for  the  generation  of  a  forecast  model 
and  hence  for  simulating  future  scenerios  of  lake  levels. 

Two  models  for  the  lake  level  simulation  were  used.  The  first  model  discussed  is 
that  of  a  deseasonahzing  filter  as  given  by  Eq.  4  followed  by  the  BURG(1967)  Max¬ 
imum  Entropy  Method  (MEM  algorithm  utilizing  the  Akaike  Information  Criterion 
(AIC)  proposed  by  AK  AIKE(  1974)  for  model  order  fitting  of  the  final  model  for  each 
lake  series.  In  the  case  of  the  BURG(1967)  MEM  algorithm,  the  AIC  criterion  is 
equivalent  to  AI<AIKE(1972)  Final  Prediction  Error  (FPE)  and  is  therefore  noted 
as  such.  The  reason  for  utilization  of  the  BURG(1967)  MEM  algorithm  is  based 
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on  the  preanalysis  of  the  spectral  shapes  of  the  deseason alized  lake  level  spectra 
which  were  investigated  both  here  and  previously  by  LIU ( 1 970)  and  are  suggestive 
of  autoregressive  (all  pole)  spectra.  The  low  frequency  content  of  the  signals  also 
suggested  that  low  order  models  might  be  reasonable  for  all  series.  This  fact  agreed 
with  the  FPE  criterion  which  had  a  minimum  final  prediction  error  for  model  order 
two  in  all  three  lake  level  series  investigated.  An  example  of  the  FPE  criterion  error 
versus  the  model  order  is  given  in  Figure  7  which  shows  a  considerable  reduction 
in  error  for  models  beyond  an  AR(1)  but  implies  that  the  error  reaches  a  broad 
plateau  beyond  a  second  order  model  suggesting  that  any  low  order  model  (two  or 
more)  might  be  used.  In  view  of  the  desirable  aspect  of  parsimony,  the  minimal  FPE 
model  order  two  was  used  for  this  model.  The  final  filter  provided  by  this  model  is 
a  combination  of  Eq.  4  followed  by  Eq.  3  with  M(B)=1.0.  Residual  testing  of  this 
model  will  be  investigated  below. 

A  second  model  investigated  was  of  the  seasonal  multiplicative  type  as  provided 
via  a  filter  consisting  of  a  seasonal  filter  of  the  type  in  Eq.  2  with  M,{B ’)  =  1.0 
and  one  seasonal  autoregressive  parameter,  followed  by  a  stationary  filter  of  Eq.  3 
form  with  M(B)=1.0  a"d  two  normal  autoregressive  parameters.  In  this  case  the 
Eq.  3  filter  was  the  BURG(1967)  MEM  algorithm  as  in  Model  1.  Residuals  from 
this  model  were  also  tested  and  will  be  discussed  below. 

RESULTS 

A  meaningful  test  of  a  model  is  provided  by  comparing  a  portion  of  the  time  history 
of  the  model  with  a  predictive  forecast  during  a  future  period  of  the  time  series  and 
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comparing  model  forecast  with  what  really  happened  during  the  same  period  of  time. 
This  approach  to  investigating  model  validity  was  used  for  the  three  lake  series  with 
the  historical  period  of  record  from  1921-1964.  The  period  of  forecast  used  was 
from  1965-1967.  It  is  noted  from  Figures  1-3  that  the  year  1964  occurred  during  a 
period  of  prolonged  low  water  level  on  the  lake  and  appeared  to  be  a  turning  point 
after  which  lake  levels  began  to  get  higher.  Forecast  series  for  the  three  lakes  for 
Model  1,  the  standardized  data  model,  based  on  the  master  gage  level  data  for  the 
lakes  are  given  in  Figures  8-10.  A  forecast  series  for  the  same  time  period  utilizing 
Model  2,  a  seasonal  multiplicative  model,  is  given  in  Figure  11.  As  can  be  seen  in 
Figures  8-11,  both  models  show  rising  water  levels  as  noted  in  nature  and  appear  to 
track  the  historical  record  within  reason  (i.e.  realizing  that  the  noise  in  the  series 
is  filtered  out).  Model  1  appears  to  provide  a  more  reasonable  approach  to  fitting 
of  the  seasonal  cycle  as  seen  in  the  historical  record,  although  the  overall  level  of 
the  series  seems  to  be  followed  more  closely  by  Model  2  (shown  for  the  case  of  gage 
5014).  The  seasonal  cycle  appears  so  prevalent  within  the  data  that  subjectively, 
at  least,  the  standardized  seasonal  reduction  model  appears  to  model  nature  more 
realistically. 

Residuals  from  both  models  were  tested  for  randomness  and  for  probability  distri¬ 
bution  of  residual  for  future  simulation  purposes.  The  original  series  were  tested  for 
normality  and  found  to  be  reasonably  normal  via  previous  work  (LIU,  1970),  and 
within  the  present  testing  framework  by  a  Kolmogrov-  Smirnov  test  (STEPHENS 
(1974)).  The  hypothesis  of  normality  was  not  rejected  with  the  Kolmogrov-Smirnov 
test  for  any  of  the  lake  level  series  at  a  five  percent  significance  level.  Based  on  this 
result,  a  normal  distribution  was  expected  for  the  residuals  since  a  Gaussian  process 
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acted  upon  by  a  linear  filter  returns  a  Gaussian  process.  This  fact  was  verified  by  a 
check  of  the  residuals  from  the  two  models.  Residuals  from  both  models  are  not  sig¬ 
nificantly  different  from  Gaussian  noise  as  determined  via  the  Kolmogrov-Smirnov 
test  at  a  five  percent  significance  level.  An  example  of  the  residual  fit  to  a  Gaussian 
distribution  is  provided  in  Figure  12  where  the  residuals  from  Model  1  with  a  two 
parameter  fit  to  gage  5014  water  levels  is  plotted  on  normal  probability  paper. 

For  simulation  purposes  (and  forecasting  purposes),  it  is  most  desirable  that  the 
model  provide  a  transformation  of  the  correlated  series  to  white  noise.  A  final  check 
on  any  model  should  investigate  whether  the  residuals  in  fact  comply  with  a  white 
noise  assumption  typically  assumed  for  simulation  purposes.  The  autocorrelation 
for  the  residuals  from  model  1  for  the  second  order  autoregressive  (MEM)  model  is 
provided  in  Figure  13  along  with  confidence  limits  for  an  assumption  of  white  noise 
at  a  five  percent  significance  level.  As  can  be  seen  in  this  figure,  the  assumption 
that  the  residuals  are  white  noise  is  well  justified.  Another  way  of  testing  this  same 
assumption  is  via  an  integrated  periodogram.  An  integrated  periodogram  of  white 
noise  follows  closely  to  a  straight  line  on  a  frequency  plot  with  the  form  7  =  af,  where 
“I”  is  a  measure  of  the  integrated  spectral  energy  and  “F  is  frequency.  The  constant 
“a”  depends  on  the  units.  Figure  14  is  an  example  of  integrated  periodogram  for 
the  residuals  from  model  1  with  a  second  order  fit  to  the  lake  level  data  of  gage 
5014.  Kolmogrov-  Smirnov  seventy-five  and  ninety-five  percent  confidence  bounds 
have  been  placed  on  the  plot  from  which  it  can  be  seen  that  the  residuals  from  the 
model  do  indeed  fulfill  the  desired  assumption  of  white  noise.  Similar  plots  for  gages 
3063  and  2030  of  the  residual  autocorrelations  and  integrated  periodograms  for  both 
models  discussed  previously  are  not  shown  here  but  provide  similar  conclusions  (i.e., 
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the  residual  series  has  been  reduced  to  white  noise). 

CONCLUSIONS 

Two  models  have  been  developed  for  use  in  simulating  future  scenerios  of  (monthly) 
average  water  levels  to  which  random  storms  can  be  superimposed  for  a  full  descrip¬ 
tion  of  the  water  level  during  an  erosion  event  or  a  future  scenerio  of  erosion  events. 
The  models  can  be  utilized  in  conjunction  with  dynamic  erosion  models  (KRAUS 
and  LARSON  (1988))  and  simulated  storm  sequences  to  determine  erosion  event 
versus  return  period.  The  models  developed  have  been  shown  to  reasonably  pre¬ 
dict  into  the  future  with  the  limited  amount  of  information  present  in  the  signal 
of  the  time  series.  The  models  have  reduced  the  correlated  seasonal  fluctuation  of 
water  level  to  a  white  noise  (uncorrelated)  series  therefore  showing  that  the  maxi¬ 
mum  amount  of  information  within  the  series  has  been  extracted  for  utilization  in 
simulating  future  scenerios  of  water  level  averages. 
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Table  1 

Parameter  Estimation  Comparison 


Period  of  Record 


1921-1954 


1955-1987 


Gage  Number 
5014 
3063 
2030 


1 .430, -.440 
1  .264, -.295 
1.399, -.431 


1 .414, -.422 
1 .149, -.167 
1 .319, -.395 


Note: 


Parameters  listed  are  the  first  and 
second  autoregressive  parameter  estimates. 


Figure  I .  Monthly  average  water  levels  -  gage  50  I  4 


Figure  5.  Auto  correlation  plot  -  gage  3063  water  levels 
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Figure  I  I .  Seasonal  multiplicative  model  forecast  -  gage  50  I  4. 
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Figure  I  4.  Integrated  periodogram  of  residuals  from  gage  50 


